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Within a class of exact time-dependent non-singular 
iV- logarithmic solutions [Mineev-Weinstein and Dawson, 
Phys. Rev. E 50, R24 (1994); Dawson and Mineev-Weinstein, 
Phys. Rev. E 57, 3063 (1998)], we have found solutions which 
describe the development and pinching off of viscous droplets 
in the Hele-Shaw cell in the absence of surface tension. 

PACS numbers 47.15. Hg, 47.20.Hw, 68.10.-m, 68.70.+w 

The process of viscous fingering in the Hele-Shaw cell, 
in which a non- viscous fluid pushes a viscous one, has at- 
tracted the attention of many physicists and mathemati- 
cians since 1958 [1]. The pressure field in the viscous 
fluid is harmonic because of D'Arcy's law, v = — Vp, 
incompressibility, and continuity, divv = 0, where v is 
fluid velocity. Because pressure is constant along the 
moving interface between the non-viscous fluid and the 
viscous fluid in the absence of surface tension, the time- 
dependent conformal technique was found to be of great 
help. The integro-differential equation 

Im(ztz^) = 1 (1) 

for the moving complex interface z(t,(f>) was derived [2]. 
Here the bar denotes complex conjugation, z t and z^ are 
partial derivatives, and the map z(t,<fi) is conformal for 
Im (j) < 0. This equation which we call the Laplacian 
growth equation (LGE) has many remarkable properties 
§- §]. In this paper, we make use of a broad class of ex- 
act time-dependent non-singular iV-logarithmic solutions 
in the channel geometry (where the interface moves be- 
tween the two parallel walls of the Hele-Shaw cell Q|) 
and in the radial geometry (where the moving contour 
is closed ||)). The equation (|I|) is valid for both of 
these geometries provided the flux of fluid is constant. 
In the channel geometry, the A^-logarithmic solution has 
the form |J] 

N 

z(t, 4>) = r(t) +^ + ^2^ log(l - Ofc(t)e-**) , (2) 

k=l 

where at = const, and |ofc| < 1. The time dependence of 
dk(t) and r(t) is given by 



{fa = z(t,i \ogcik) 
= t - logSfe + Y,i=i a i lo g(! - a fc a;) = const, (3) 
t + C ={l - ^J2fLi a l) T + | E/Il lo g( a > 

where k — 1,2, ...,7V and C are constants in time. 
The Hele-Shaw cell width is chosen to be 27r. Hence 
z(2tt) — z(0) = 2iri. This is for periodic boundary condi- 
tion, while for the no-flux boundary conditions the same 
formulas (|^-||) are valid, but the sums in (|[]|) contain 
both terms with ctk and and with uu and unless 
both cvfc and ak are real. 

The A^-logarithmic solution in radial geometry is 

z(t,d>) = r(t)e** + £ a fc log(_--l), (4) 

where X^fc=i a k = 0, cvfe = const, and |ojt| < 1. 
Time dependence of afc(i) and r(t) is given by 

{Pk = z(t,i loga fc ) 
— const 
2i + C = r 2 -r£f =1 ^. 

For a broad range of {ev^'s} (which we do not specify 
here) the solutions (||) are free of finite-time singularities 
(cusps) 1-1. 

In this paper we present a subclass of solutions (g) 
which describes the change of topology of the moving 
interface which occurs in the process of droplet forma- 
tion of the viscous fluid. The moving interface contin- 
ues to evolve in accordance with the equation (|l[). We 
have found the analytic continuation from the domain 
immediately before the topological breakdown occurred 
to the domain immediately after the droplet formation. 
This continuation allows us to study the interface dy- 
namics which includes an arbitrary number of topological 
changes of this kind in the absence of surface tension. We 
believe that inclusion of surface tension in these studies 
will shed light on which topological changes are physi- 
cally allowed and which are forbidden. The goal of this 
article is to report a new class of solutions which describes 
droplet formation in the absence of surface tension. 

Let us consider the following 3-logarithmic solution of 
the LGE (jl|) in channel geometry: 
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z(t, (j>) = r(t) + i<fi + a log ■ 



ai4> 



a(t) 



,i<p 



-b(t) 



(5) 



where a > 0, < 6(0) < a(0) < 1, and the third logarith- 
mic singularity, corresponding to the term i<j>, is zero. It 
follows from that the time evolution of ait) and b(t) 
obeys the equations 



/3i = £ - f log f + a log - log a 
/3 2 = t - | log f + a log ~ log6 



(0) 



Both a(f) and b(t) reach the unit circle simultaneously at 
the moment, t* 1 which is easy to calculate: 

Since a and b are close to 1, then deviations of a and 
b from 1, defined as p = 1 — a and q = 1 — b, are very 
small, that is < p = 1 — a < q = 1 — b « 1. Thus the 
last system in the leading order is 



Pi = t* + a log c — a log 
f3 2 = t* + a log c + a log 



. c+c 



(7) 



where c 2 = lim{_>t» p/g is a finite constant (assuming 
that t* is finite). The solution of (0) is 



c = A - ^A 2 - 1 

t* = 0x - a log(l - v 7 ! - e^i-Z^A*, 



(8) 
(9) 



where A = eC^-ft)/ 2 ". 

It is also easy to calculate the behavior of a(t) and b(t) 
immediately prior to the moment t*: 

Defining St — t* — t and expanding a(t) and b(t) near 

t* 



■ p St - Pl (5t) 2 
q 5t - qi(St) 2 , 



(clearly then, that po/qo = c 2 ), we rewrite 
linear terms with respect to St: 



= a(— — Pl+Po 
^po qi+qo 

V oi +an an I ^ u v 2 



„ (a qa-pa 

Pol ^ 



■ 91+90 90 

Keeping in mind that c 2 
this last system that 

J Po 



qa+Po 

qo-po 

2 90+Po 



1) 
1) 



(10) 
keeping 

(11) 



Po/qo (see above) we find from 



(l+a/2)+c a (l-Q/2) 
9° = (l+a/2)+c a (l-a/2) ' 



(12) 



where c is given by 
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FIG. 1. The front, moving to the right, forms a viscous 

droplet, and then leaves it behind. The numbers 1-7 denote 

7 distinct snapshots of the advancing front. 

It is interesting to see how the advancing front moves 
in this case. Figure [l] shows that the initially almost 
planar interface (6(0) << 1) advances until a groove of 
width na is created at the stagnation point with com- 
plex coordinate (3\ — a log 2. Then later this groove be- 
comes thinner and thinner until at the moment t* it dis- 
appears completely. Since both a(t) and b(t) reach the 
unit circle simultaneously, (||) contains a singular term 
log[(e I * — l)/(e 1( ^ — 1)] which has a removable singularity 
at 1. This term equals zero exactly at this moment of 
time. If we formally consider t > t* in (6) we see that 
both a(t) and b(t) continue their evolution outside the 
unit circle. This is clearly an unphysical part of the solu- 
tion (||) , which violates the framework of validity of this 
solution, which demands that the exterior of the unit cir- 
cle should be free of singularities except at infinity. We 
also note that the range of the variable <fi parametrizing 
the moving curve, corresponding to the groove of viscous 
liquid which is left behind the advancing front, becomes 
negligibly small when the "neck" of this groove becomes 
very thin. This allows us to make the analytical con- 
tinuation from the moment t*_ immediately prior to the 
disappearance of the neck and pinching off the viscous 
droplet behind the moving interface to the moment t* + 
that is immediately after this event occurs. Once both 
singularities a(t) and b(t) reach the unit circle at t = t* , 
the last term in (0) vanishes and the subsequent evolu- 
tion ignores this term. The interface in this case becomes 
a planar front described by the equation 



z(t,ct>)=T(t)+i<p, 



(13) 



where r = t in this case. In the physical plane, this an- 
alytical continuation at the moment t — t* shows that 
the viscous droplet, which was formed from the original 
groove at the stagnation point j3\ — a log 2, pinches off 
from the advancing interface and is left behind the mov- 
ing front. The moving front for t > t* is described by the 
same solution but without the pair of logarithms. This 
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process can be seen in Figure [j]. The shape of the droplet 
does not change after t — t* . This is because of two lim- 
itations of the present theory: both the viscosity of the 
less viscous liquid and the surface tension were chosen to 
be zero. This leads the pressure along the droplet bound- 
ary to be a constant, and then, by virtue of the maximum 
theorem for harmonic functions, the pressure should be 
constant throughout the interior of the droplet. The pres- 
sure gradient acting on the droplet is zero. Consequently 
the droplet is at rest in the non-viscous domain. In the 
real situation, where neither the second liquid viscosity 
nor the surface tension is zero, the viscous droplet will 
be dragged by the flow and its shape will continue to be 
modified. 

The shape of the droplet described by (H) has a simple 
form: this is an oval with a horizontal length of the order 
of 02 — /3i + 2odog2 and the vertical size about ira. The 
droplet's shape can be, however, much more complex. It 
can be described by the N-logarithmic solution (jsj) such 
that the ctk corresponding to the minimal \au\ is cho- 
sen to be — ai, where the prime in the sum indicates 
summation over all I from 1 to N, but omitting I = k. 

In Q the geometrical interpretation of the constants 
of motion and {Pk} have been found. Namely, 

(5k — a,k log 2 is the complex coordinate of the stagna- 
tion point near which the k th groove with parallel walls 
originates. And ir\ak\ and arg(afc) are the width and 
the angle with respect to the horizontal axes respec- 
tively. This geometrical interpretation appeared to be 
in the excellent agreement with known real and numeri- 
cal experiments. But this interpretation has a constraint: 
arg(afe)| < n/2. In this work we consider ol^ which vi- 
olates the last inequality. For such ak the /3* — afclog2 
is the complex coordinate of the point where the pre- 
viously formed groove starts to shrink. Clearly, for the 
simplest example with two logarithms, where a± > and 
a2 < 0, the width of the groove which was originally noti 
becomes after shrinking n(ai +02). Of course, this con- 
sideration is valid only when otx + ct2 > 0. In the case 
of equality the width vanishes completely, it corresponds 
to the formation of the droplet described above. We do 
not understand on this stage why and when the terms 
with Im(afc) < are physically realizable. We hope that 
future studies will elucidate this point. 

Let us derive the droplet's neck dynamics w(t) prior to 
pinch off in the example d) discussed above. The width 
w(t) of the neck is defined as 

w(t) = 2y(t,^*(t)) = 2Imz(t,0* (0, (14) 

where (f>* is a non-zero solution of the equation 

dy/dcj, = 0, (15) 

as it is clear fro m da ) and from Figure [l]. Differentiating 
(||) and solving (|15|), having in mind that both a and b 
are close to 1, we obtain </>* = y/ a(a — b) in the leading 



order with respect to (a — b) which is small. Then from 
(pi) and (|) it follows that 



w(t) = 2[(fr* + a(arctan 



cos 4>* — a 



arctan ■ 



COS ( 



(16) 



Replacing here 
obtain 



by its value found above we finally 



w{t) = A^a{a -b) = Va(l - c 2 )q (t* - t), (17) 

where constants c, qo, and t* are given above. 

Here we considered only the simplest configuration d) 
leading to the formation of a single droplet. Clearly, 
there is a broad subclass of solutions (|J) and (|4|) which 
describes in a similar fashion formation and pinching 
off of an arbitrary number of viscous droplets with dif- 
ferent shapes. The same analytic continuation, per- 
formed above for the simple case N = 3, holds for 
this generalized case. Again, this continuation contains 
the elimination of the removable singularity by equating 
(e 1 ^ — l)/(e 1 ^ — 1) to 1 at the moment of formation of 
each droplet, thus decreasing the total number of loga- 
rithms in the description of the moving interface, as done 
above for N — 3. All these considerations hold for both 
the channel geometry and for closed interfaces. 

Previous studies on topological changes in the viscous 
flow in the Hele-Shaw channel || considered surface ten- 
sion as a significant factor. Because surface tension is 
zero in the present paper and because we believe that 
there are many different regimes of pinching off within 
classes of solutions (||) and ([|), we do not compare our 
findings with results in ||. There is also a paper on the 
change of topology in the absence of surface tension [Q, 
which states that a change of topology is inevitable in 
the interior (suction) problem if the sink is not located 
at the center of mass of the viscous domain. Here we 
may say that we know how to describe possible topolog- 
ical changes in detail. What we do not know is which 
of the solutions obtained above are physically realizable 
and which are not. We also do not know at this stage 
what will be the corrections to the dynamics of droplet 
formations caused by non-zero surface tension. We hope 
that future studies will answer these questions. 
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